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For a theory with a pseudo scalar coupling 4>FF and in the case that there 



J> ■ is a constant electric or magnetic strength expectation value, we compute the 

o . 

interaction potential within the structure of the gauge-invariant but path- 
. dependent variables formalism. While in the case of a constant electric field 

strength expectation value the static potential remains Coulombic, in the 
case of a constant magnetic field strength the potential energy is the sum of a 
Yukawa and a linear potentials, leading to the confinement of static charges. 



PACS number(s): ll.10.Ef, 11.15.-q 

I. INTRODUCTION 

One of the key issues facing QCD is understanding confinement of quarks and gluons. In 
fact, a linearly increasing quark- ant iquark pair static potential provides the simplest criterion 
for confinement, although unfortunately there is up to now no known way to analytically 
derive the confining potential from first principles. However, as is well known, it has been 
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approached from many different techniques and ideas, like lattice gauge theories [1], non 
perturbative solutions of Schwinger-Dyson's equations [2]. Other authors also associate 
confinement with the existence of a non trivial vacuum structure, where the chromomagnetic 
field strength acquires a non zero expectation value [3]. 

In fact non vanishing expectation values can have dramatic consequences in everything 
that concerns the infrared properties of a theory. An interesting model where these effects 
have been studied to a certain extent is the " axion-gauge field" system, where a scalar field 
(j) (the "axion") is coupled to gauge fields via the interaction term 

Cj = fyeT^F^F^. (1) 

This theory experiences mass generation when <p develops a space dependent expectation 
value [4] and tachyonic mass generation when a time dependent expectation value appears [5]. 
Mass generation is also achieved when the gauge field takes a magnetic type expectation 
value [6]. If F^ v takes an electric type expectation value, tachyonic mass generation takes 
place [6]. Thus, in order to gain further insight into the physics presented by this theory, 
in this paper we will focus attention on the static potential between charged fields. The 
purpose here is to investigate the effects of the external expectation value field strength on 
the interaction energy. 

The interaction energy between static charges is a tool of considerable interest which is 
expected to provide the foundation for understanding confinement, and its physical content 
can be understood when a correct separation of the physical degrees of freedom is made. 
Previously, we proposed a general framework for studying the confining and screening nature 
of the static potential in gauge theories in terms of the gauge- invariant but path-dependent 
field variables [7] . An important feature of this methodology is that it provides a physically- 
based alternative to the usual Wilson loop approximation. When we compute in this way 
the static potential for the model described in [6], which contains the term (1), in the 
presence of an external field strength which can be either electric or magnetic, the result of 
this calculation is rather unexpected in the magnetic case: It is shown that the interaction 
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energy is the superposition of a Yukawa and a linear potentials, that is, the confinement 
between static charges is obtained. On the other hand, in the case of a constant electric field 
strength expectation value the static potential remains Coulombic, that is, the interaction 
energy does not exhibit any sensitive modification. Actually, the linear confining potential 
seems to be associated only with the magnetic field strength expectation value. 

II. INTERACTION ENERGY 

Before going to the derivation of the interaction energy, we will describe very briefly the 
model under consideration. We start from the following effective Lagrangian [6]: 

£ = -\f^r - Y^ uaP (^) £PaiS <^> teuhrfj'*' (2) 

where (F^) represents the constant classical background (which is a solution of the classical 
equations of motion), and uia is the mass for the axion field. Here, f^ u = d^A u — d u A fM 
describes a small fluctuation around the background. We also mention that the above 
Lagrangian arose after using e^ al3 (F^) (F^) = (which holds for a pure electric or a pure 
magnetic background), and integrating out the axion fields 

By introducing e^ 13 (F^ v ) = v a/3 and e pal& (F pr7 ) = v 1& , it follows that the expression 
(2) can be rewritten as 

1 a 2 1 

C = — LJ^ - —v a ?f a0 T^f-ts, (3) 

4 ^ 16 1 a + m\ 7 v ; 

still, the tensor v af3 is not arbitrary, but must satisfy e^^v^Vafs = 0. 

A. Magnetic case 

As stated, our main objective is to calculate the interaction energy in the v° l ^ and 
v lj = case (referred to as the magnetic one in what follows), following the conventional 
path via the expectation value of the Hamiltonian in the physical state |<E>), which we will 
denote by (H)$. The Lagrangian (3) then becomes 
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£ = --J^r - ^v 0i foi^-^v° k fok - A)A (4) 
4 16 □ + m A 

where J° is an external current, (//, v — 0, 1, 2, 3) and (i, A: = 1,2,3). 

We now proceed to obtain the Hamiltonian. For this we consider the Hamiltonian for- 
mulation of this theory. The canonical momenta obtained from (4) are 

n° = o, (5) 

and 

n* = DijEj, (6) 



where Ej_ = F i0 and = ySij — Y v io n+m 1 v i°J • Since D is a nonsingular matrix (det £) = 
1 ~~ ^5^5- 7^ 0) with v 2 = v l0 v l °, there exists the inverse of D and from Eq.(6) we obtain 

E t = — — (<y i7 - det D + — ^rv 7 | n.-. (7) 

detD \ lJ 8 D + m 2 A 3 j J K ' 



The canonical Hamiltonian corresponding to (4) is 



h c = J d 3 x j-A) (^ir - J°) + in 2 + 



s 2 (v-n) 2 , i. 



+ o B ; (8) 



16(D + M 2 ) 2 

where M 2 = m 2 A — ^-v 2 an d B is the magnetic field. Demanding that the primary con- 
straint Il = be preserved in the course of time, one obtains the secondary Gauss law 
constraint of the theory as Ti (x) = dill 1 — J° = 0. The preservation of r\ for all times 
does not give rise to any further constraints. The theory is thus seen to possess only 
two constraints, which are first class, therefore the theory described by (4) is a gauge- 
invariant one. The extended Hamiltonian that generates translations in time then reads 
H = He + / d 3 x (co (x) n (x) + c\ (x) T± (x)), where Co (x) and c\ (x) are the Lagrange mul- 
tiplier fields. Moreover, it is straightforward to see that A (x) = [A (x) , H] — cq (x), which 
is an arbitrary function. Since n° = always, neither A nor n° are of interest in describing 
the system and may be discarded from the theory. Thus the Hamiltonian takes the form 

H= [ d 3 x[-U 2 + 9 - , (V - n) l + -B 2 + c(x) (d t W -J )), (9) 
J [2 16(D + M 2 ) 2 v ; V // ' K ' 



where c(x) = c\(x) — A (x). 

To quantize the theory using Dirac's procedure [8] we introduce a supplementary condi- 
tion on the vector potential such that the full set of constraints becomes second class. For 
this purpose, we could choose, for example, the gauge-fixing condition [7] 

1 

T 2 (x) = J dz u A u (z) = J d\x l Ai (Xx) = 0, (10) 

c (x o 

where A (0 < A < 1) is the parameter describing the spacelike straight path x l — £ l + 
A (x — and £ is a fixed point (reference point). There is no essential loss of generality 
if we restrict our considerations to = 0. In this case, the only nonvanishing equal-time 
Dirac bracket is 

i 

{A t (x) , IF (y)}* = 6*6® (x-y)-dfj dXx>5^ (Xx - y) . (11) 

o 

In passing we recall that the transition to quantum theory is made by the replacement of 
the Dirac brackets by the operator commutation relations according to 

{A B}* -> (-i) [A, B] . (12) 

We are now in a position to evaluate the interaction energy between pointlike sources in 
the model under consideration, where a fermion is localized at yf and an antifermion at y. 
From our above discussion, we see that (if)$ reads 

(H)^ = ($| [ d'xllu 2 + t- ^ ' Si + i B 2 | |$) . (13) 



2 16(D + M 2 ) 2 
Next, as was first established by Dirac [9], the physical state can be written as 

|$) = |* (y) * (yf)) = i> (y) exp J dz l A t (z)j ^ (yf) |0) , (14) 

where |0) is the physical vacuum state and the line integral appearing in the above expression 
is along a spacelike path starting at yf and ending at y, on a fixed time slice. From this we 
see that the fermion fields are now dressed by a cloud of gauge fields. As mentioned before, 
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the fermions are taken to be infinitely massive (static). Consequently, we can write Eq.(13) 

as 

m .. m j*$?-£g*} m , (15) 

with = -V 2 . 

From the foregoing Hamiltonian discussion, we first note that 

n,(x) * (y) * (y')) = * (y) * (y') n t (x) |0> + e F 4*8™ (z - x) |*> . (16) 

Combining Eqs.(15) and (16), we have 

{H) 9 = (H) + V 1 + V 2 , (17) 

where (H) = (0\H\0). 
The V\ term is given by 

V, = e -\ y dz'df [ y dz%G (z',z), (18) 

A Jy Jy 

where G is the Green function 

1 -M|z'-z| 
47T |Z — Z| 

By means of Eq.(19) and remembering that the integrals over z % and z\ are zero except on the 
contour of integration, the term (18) reduces to the Yukawa-type potential after subtracting 
the self-energy terms, that is, 

2 -M\y-y'\ 

Vl = ~~ jr. (20) 

4tt |y-y'| 

We now come to the V 2 term, which is given by 

V 2 = F dz H F dz'Giz', z) . (21) 

2 Jy Jy 

In order to compute V 2 , we make use of the Green function (19) in momentum space 

1 e -M|z'-z| d 3 k ik-(»'-«) 

(22) 



An \z'-z\ J (27r) 3 k 2 + M 2 ' 
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Thus, by employing relation (22) we can reduce Eq.(21) to 

r d 3 k 1 1 

V 2 = e 2 m 2 A o [1 - cos (k • r)] — (23) 

A J (2tv) 3 1 V n (k 2 + M 2 ) (n-k) 2 V ; 

where n = | y Iy/| is a unit vector and r = y — y' is the relative vector between the quark 

and antiquark. Since n and r are parallel, we get accordingly 

where £>r denotes the momentum component perpendicular to r. We may further simplify 
Eq.(24) by doing the kr integral, which leads immediately to the result 

v *-^Jw ll -«*^H 1 + *Tm)' (25) 

— OO \ I / 

where A is an ultraviolet cutoff. We also observe at this stage that similar integral was 
obtained independently in Ref. [10] in the context of the dual Ginzburg-Landau theory by 
an entirely different approach. 

Now, we move on to compute the integral (25). To this end it is advantageous to 
introduce a new auxiliary parameter e by making in the denominator of the integral (25) 
the substitution k 2 — > k 2 + e 2 . This allows us to obtain a form more comfortable to handle 
the integral. Hence we evaluate lim e ^ V 2 , that is, 

* - S3 * - S3 ^ 7 wh) I 1 - cos <^>1 ( l + V^w) ■ (26) 
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The integration on the A; r -complex plane yields 
Taking the limit e — > 0, expression (27) then becomes 

^^! ly _ y11n ( 1 + ^). (28) 

This, together with Eq.(20), yields finally 



where L = |y — y'|. 

It is worth noting here that this is exactly the result obtained in Ref. [10] in the context of 
the dual Landau-Ginzburg thaeory. But we do not think that the agreement is an accidental 
coincidence. Also, the massive Abelian antisymmetric tensor gauge theory displays the same 
behavior [11,13]. In other words, there is a class of models which can predict this interaction 
energy. 

B. Electric case 

We now want to extend what we have done to the case v° l = and v lj ^ (referred to 
as the electric one in what follows). In such a case the Lagrangian reads 

C = -\fuuF u - —v^h l -^v kl f kl - A J°, (30) 

(yu, v — 0, 1, 2, 3) and (i, j, k, I — 1, 2, 3). 

The above Lagrangian will be the starting point of the Dirac constrained analysis. The 
canonical momenta following from Eq.(30) are IP = which results in the usual primary 
constraint II = and II* = f l °. Defining the electric and magnetic fields by E % = F t0 and 
B l = j£^ k Fjk, respectively, the canonical Hamiltonian assumes the form 

H C = J d 3 x | ^E 2 + l -B 2 + ^s ijm s kln v^B m ^-^v kl B n - A (dJV - J°) j . (31) 

Time conservation of the primary constraint leads to the secondary constraint Ti(x) = 
<9jIP — J° = 0, and the time stability of the secondary constraint does not induce more 
constraints, which are first class. It should be noted that the constrained structure for the 
gauge field is identical to the usual Maxwell theory. Notwithstanding, in order to put our 
discussion into context it is useful to summarize the relevant aspects of the analysis described 
previously [12]. In view of this situation, we pass now to the calculation of the interaction 
energy. 

Following our earlier procedure, we will compute the expectation value of the Hamiltonian 
in the physical state |<&) (Eq. (14)). That is, 
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<#> $ H$l/^{^E 2 }|$>. (32) 
Taking into account the above Hamiltonian structure, the interaction takes the form 

(H)* = (H) + € - r dz* f y dz[ 5^(z - z'), (33) 

Z Jy Jy 

where {H) — (0 | H | 0). Once again, the integrals over zi and z\ are zero except on the 
contour of integrations, accordingly one obtains the following interaction energy: 

v = y*|y-y'l , (34) 

where k = 5^(0). This expression shows that special care has to be exercised in order to 
clarify the appearance of this peculiar result, as was discussed elaborately in [12]. It may 
be recalled, however, that the origin of the divergence is quite clear, so that it is possible 
to extract the Coulomb potential from the infinite contribution. Notice that the origin of 
the divergent factor k is due to the fact that the thickness of the string is nonvanishing 
only on the contour of integration. We recall that a suitable examination of the term 
Y J d 3 x (jy ' dzid^ (x — z)) reproduces exactly the expected Coulomb interaction between 
charges after subtracting the self-energy term [12], hence Eq.(34) reduces to 

V(L) = -£i (35) 

where L = | y — y ' | . 

III. FINAL REMARKS 

We briefly summarize the results obtained so far. By using the gauge-invariant but 
path-dependent formalism, we have studied the static potential for the system consisting 
of a gauge field interacting with a massive axion field in the case when there are nontrivial 
constant expectation values for the gauge field strength F^ u . The constant gauge field 
configuration is a solution of the classical equations of motion. 
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While in the case when {F^J) is electric-like no unexpected features are found, we find that 
the case when (F^) is magnetic-like is totally different. In fact, when (F^) is magnetic- 
like, the potential between static charges displays a Yukawa piece plus a linear confining 
piece. Unexpectedly, a confining potential between static charges appears in this case. It 
is interesting to note that the requirement that (F^ u ) be magnetic in order to get confining 
behavior coincides with ideas concerning the nature of the QCD vacuum, where a nontrivial 
magnetic field strength must be present in the vacuum [3]. 

In the non-Abelian generalization of this model similar effects should appear. 

Further investigations of the relations between our work and the magnetic models of the 
QCD vacuum in Ref. [3] have to be performed. Also one should address the question of 
what is the physical origin of the axion field, it is probably some kind of a bound state, if 
our theory represents an effective approach to QCD. 
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